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1. Introduction 

Let (M, g) be a Riemannian manifold of dimension n > 2. Strichartz estimates 
are a family of space time integrability estimates on solutions u(t, x) : (— T, T) x 
M ^ C to the wave equation 

(1.1) dtu{t,x) - Agu{t,x) ^0, u{0,x) = f{x), dtu{0,x) ^ g{x) 

where Ag denotes the Laplace-Beltrami operator on {M,g). Local homogeneous 
Strichartz estimates state that 

(1-2) \\u\\lp{{~T,T);Li{M)) < C {\\f\\Hi(M) + ll.9!l_ff^-i(j\/)) 

where denotes the Sobolev space over M of order 7, and 2 < p < 00, 
2 < g < 00 satisfy 

Inn 2 n — I n — 1 

1.3 _ + _= _ + < 

p q 2 P q 2 

Estimates involving q = 00 hold when (n,p, q) 7^ (3, 2, 00), but typically require the 
use of Besov spaces. 

Strichartz estimates are well established on flat Euclidean space, where M — R" 
and gij — Sij. In that case, one can obtain a global estimate with T = 00; see for 
example Strichartz [77], Ginibre and Velo [S], Lindblad and Sogge [Tfr, Keel and 
Tao [14, . and references therein. However, for general manifolds phenomena such as 
trapped geodesies and finiteness of volume can preclude the development of global 
estimates, leading us to consider local in time estimates. 

If M is a compact manifold without boundary, finite speed of propagation shows 
that it suffices to work in coordinate charts, and to establish local Strichartz esti- 
mates for variable coefficient wave operators on R". Such inequalities were devel- 
oped for operators with smooth coefficients by Kapitanski and Mockenhaupt- 
Seeger-Sogge |18) . In this context one has the Lax parametrix construction, which 
yields the appropriate dispersive estimates. Strichartz estimates for operators with 
C^'^ coefficients were shown by the second author in 21j and by Tataru in '29\, 
the latter work establishing the full range of local estimates. Here the issue is more 
intricate as the lack of smoothness prevents the use of the Fourier integral operator 
machinery. Instead, wave packets or coherent state methods are used to construct 
parametrices for the wave operator. 
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In this work, we consider the estabhshment of Strichartz estimates on a manifold 
with boundary, assuming that the solution satisfies either Dirichlet or Neumann 
homogeneous boundary conditions. Strichartz estimates for certain values of p, q 
were established by Burq-Lebeau-Planchon [5j using results from 25J; our work 
expands the range of indices p and q, and includes new estimates of particular 
interest for the critical nonlinear wave equation in dimensions 3 and 4. Our main 
result concerning Strichartz estimates is the following. 

Theorem 1.1. Let M he a compact Riemannian manifold with boundary. Suppose 
2<p<oo, 2 < q < oo and {p, g, 7) is a triple satisfying 

Then we have the following estimates for solutions u to (jl.ip satisfying either 
Dirichlet or Neumann homogeneous boundary conditions 

(1-5) \\u\\lp([-T.T\;Li{M)) < C ( II/I|h^(A/) + \\g\\Hi-^(M)) 

with C some constant depending on M and T . 

A lemma of Christ-Kiselev [7] allows one to deduce inhomogeneous Strichartz 
estimates from the homogeneous estimates. In the following corollary, (r',s') are 
the Holder dual exponents to (r, s) , and the assumptions imply that a homogeneous 
{H^-'^,H-^) UL" holds. 

Corollary 1.2. Let M be a compact Riemannian manifold with boundary. Suppose 
that the triples {p, g, 7) and (r', s' , 1— 7) satisfy the conditions of Theorem \l.l[ Then 
we have the following estimates for solutions u to (jl.ip satisfying either Dirichlet 
or Neumann homogeneous boundary conditions 

\W\\li'{[-T,T]-Li{M)) < C ( \\f\\Hi{M) + I|5||h-'-1(M) + \\P\\l^{[-T,T\-L'(M))) 

with C some constant depending on M and T . 

For details on the proof of CoroUarv l 1 . 21 using Theorem ll.ll and the Christ-Kiselev 
lemma we refer to Theorem 3.2 of [24 , which applies equally well to Neumann 
conditions. 

By finite speed of propagation, our results also apply to noncompact manifolds, 
provided that there is uniform control over the size of the metric and its derivatives 
in appropriate coordinate charts. In particular, we obtain local in time Strichartz 
estimates for the exterior in K" of a compact set with smooth boundary, for metrics 
g which agree with the Euclidean metric outside a compact set. In this case one 
can obtain global in time Strichartz estimates under a nontrapping assumption. 
We refer to [H] for the case of odd dimensions, and Burq [3] and Metcalfe [T7] for 
the case of even dimensions. See also [TI] . 

For a manifold with strictly geodesically-concave boundary, the Melrose- Taylor 
parametrix yields the Strichartz estimates, for the larger range of exponents in 
(1.3) (not including endpoints) as was shown in [23j . If the concavity assumption 
is removed, however, the presence of multiply reflecting geodesies and their limits, 
gliding rays, prevent the construction of a similar parametrix. 
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Recently, Ivanovici [12] has shown that, when n — 2, (1.5) cannot hold for the 
full range of exponents in (1.3). Specifically, she showed that if M C is a 
compact convex domain with smooth boundary then (1.5) cannot hold when g > 4 
if 2/p + 1/q = 1/2. It would be very interesting to determine the sharp range of 
exponents for (1.5) in any dimension n > 2. 

The Strichartz estimates of Tataru [2^ for Lipschitz metrics yield estimates in 
the boundary case, but with a strictly larger value of 7. The approach of p5] 
involves the construction of parametrices which apply over short time intervals 
whose size depends on frequency. Taking the sum over such sets generates a loss of 
derivatives in the inequality. 

These ideas influenced the development of the spectral cluster estimates for man- 
ifolds with boundary appearing in [25j . Such estimates were established through 
squarefunction inequalities for the wave equation, which control the norm of u{t, x) 
in the space L'^{M; L^{—T,T)). These spectral cluster estimates were used in the 
work of Burq-Lebeau-Planchon [5] to establish Strichartz estimates for a certain 
range of triples {p, q, 7). The range of triples that can be obtained in this manner, 
however, is restricted by the allowed range of q for the squarefunction estimate. 
In dimension 3, for example, this restricts the indices to p,q > 5. In [5] simi- 
lar estimates involving W'"^ spaces were also established, and used in conjunction 
with the Strichartz estimates and boundary trace arguments to establish global 
well-posedness for the critical semilinear wave equation for n — 3. In the last two 
sections of this paper we shall present some new results concerning critical semi- 
linear wave equations. Specifically, we shall obtain local well-posedness and global 
existence for small data when n = 4, as well as a natural scattering result for n — 3. 

The approach of this paper instead adapts the proof of the squarefunction in- 
equalities in [25j . We utilize the parametrix construction of that paper, and estab- 
lish the appropriate time-dispersion bounds on the associated kernel. This allows 
us to obtain the Strichartz estimates for a wider range of triples, including, for 
example, the important L'^((—T,T);L^'^{M)) estimate in dimension 3, and the 
L^{{-T,T);L^{M)) estimate in dimension 4. 

The key observation in [25] is that u satisfies better estimates if it is microlocal- 
ized away from directions tangent to dM than if it is niicrolocalized to directions 
nearly tangent to dM. This is due to the fact that one can construct parametrices 
over larger time intervals as one moves to directions further away from tangent to 
dAI. More precisely, the parametrix for directions at angle « 6 away from tan- 
gent to dM applies for a time interval of size 9, which would normally yield a 
6'-dependent loss in the estimate. However, this loss can be countered by the fact 
that such directions live in a small volume cone in frequency space. For sub-critical 
estimates, i.e. where strict inequality holds in the second condition in (|1.3p . this 
frequency localization leads to a gain for small 6. The restriction on p and q in 
Theorem 1 1 . 1 1 arises from requiring this gain to counteract the loss from adding over 
the disjoint time intervals on which one has estimates. Hence, while the range 
of p and q in our theorem is not known to be optimal, the restrictions are naturally 
imposed by the local nature of the parametrix construction in [25j . 
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Notation. The expression X < Y means that X < CY for some C depending 
only on the manifold, metric, and possibly the triple (p, q, 7) under consideration. 
Also, we abbreviate LP{r,Li{U)) by LPL'^^I x U). 

2. Homogeneous Strichartz Estimates 

The proof of Strichartz estimates is a direct adaptation of the proof of square- 
function estimates in [25j . The difference is that Strichartz estimates result from 
time decay of the wave kernel, whereas squarefunction estimates result from de- 
cay with respect to spatial separation. Consequently, in |25j the wave equation 
was conically localized in frequency so as to become hyperbolic with respect to a 
space variable labelled xi, and the equation factored so as to make xi the evolution 
parameter. 

In order to maintain the convention that xi is the evolution parameter, in this 
section we set xi — t, and will use x' = {x2, ■ ■ ■ , a;„+i) to denote spatial variables 
in M". Thus x — {xi^x') is a variable on M-'-+". 

We work in a geodesic-normal coordinate patch near dM in which a;„ > equals 
distance to the boundary (the estimates away from dM follow from [l3j and [18]). 
The coefRcients of the metric gij{x') are extended to a;„ < in an even manner, and 
the solution u{x) is extended evenly in the case of Neumann boundary conditions, 
and oddly in Xn in case of Dirichlet conditions. The extended solution then solves 
the extended wave equation on the open set obtained by reflecting the coordinate 
patch in Xn- 

Setting a^^{x) = ^detgij{x'), we now work with an equation 

n+l 

Dia'^{x')Dju{x)=0 

on an open set symmetric in a;„. A linear change of coordinates, and shrinking the 
patch if necessary, reduces to considering coefficients a*-' (x) which are pointwise 
close to the Minkowski metric on the unit ball in M^+", and defined globally so as 
to equal that metric outside the unit ball. 

Following f5F, §2], the solution u is then localized in frequency to a conic set 
where |^'| ~ |^i|. On the complement of this set the operator is elliptic, and the 
Strichartz estimates follow from elliptic regularity and Sobolev embedding. As in 
section 7 of one uses the fact that the coefficients are smooth in all variables 
but Xn, and Sobolev embedding can be accomplished using at most one derivative 
in the x„ direction. 

The next step is to take a Littlewood-Paley dyadic decomposition u — X]fc°=i 
with Mfc localized in frequency to shells |^'| ~ 2'^. One lets a^{x) denote the 
coefficients frequency localized in the x' variables to |^'| < 2*^, and factorizes 

n+l 

where p/c(x, « |^'|. Just as in [25, §2], (and the higher dimensional modifications 
in [25, §7]), Theorem ll.il is reduced to establishing, uniformly over A = 2*^, bounds 
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of the form 

(2.1) ||ua||l£^l^,(|,|<i) < Xmiuxh^L- + \\Fx\\l-) , Dmx - Px{x,D')u = Fx . 

Here, Px{x,D') — ^px{x,D') + ^px{x, D')* , and the symbol px{x,^') can be taken 
frequency locahzed in x' frequencies to |^'| < A, and px{x,£,') = \C\ if IC'I 9^ 

The setup is now the same as in |25j . and the reductions of §3-§6 of that paper, 
specifically their n-dimensional analogues of §7, apply directly. This starts with a 
decomposition ux = corresponding to a dyadic decomposition of uxiS,) in 

the ^„ variable to regions ^„ e [2-^-'^X, 2"^+^] where A-^/^ < 2-0 < 1. 

If 2^^ > i, corresponding to non-tangential reflection, then the estimates will 
follow as the case for 2~^ = i, so we restrict attention to the case 2"-' < i. Since 
l^'l A, this implies that some remaining variable is ~ A, and after rotation we 
assume that Uj{xi,£,') is supported in a set 

■ 61+1 ~ A, I < cA, j = 2, . . . , 71 - 1, and ^„ « 9j\} 

where X~^/^ < Oj < I. 

The proof establishes good bounds on the term uj over time intervals of length 
9j. Precisely, let Sj^k, \k\ < 0~^, denote the time slice Xi G [ke6j,{k+ l)e0j]. In 
analogy with pSl Theorem 3.1], we establish the bound 

(2.2) ll%IU£,L2,(S,,0^^"<''''c.,fc 

where Cj^k satisfies the nested summability condition [251 (3.1)], and where 

in-m-i)-h (--2)(i-|)<| 

h-h (--2)(i-i)>f 

Adding over the 6J^ disjoint slabs intersecting \xi \ < 1, the simple uniform bounds 
on Cj^k yield 

lkjllLS,L^,(|.,|<i) < A^ef Wuxh^L^ + WFxWl^) . 

The 9j take on dyadic values less than 1, and provided a{p, q) > 1/p, one can sum 
over j to obtain (|2.ip . In case a{p, q) ~ l/p one can also sum the series, using the 
nested summability condition |25l (3.1)], together with the branching argument on 
[25l page 118], to yield (|2.ip . Note that the restrictions on {p,q) in Theorem 11.11 
are precisely that a{p, q) > l/p. 

Estimate (|2.2p is established through the parametrix construction from [52], 
together with the use of the V2 spaces of Koch-Tataru [15 . Precisely, one rescales 
jji+n g^^^ considers the symbol 

q{x,^')^9,pj{ejx,ej'a, 

where pj is such that q{xiX,S.') is supported in |C| < cfi^^'^, where /i = 9jX is the 
frequency scale at which Uj{djx) is localized. Fix u{x) = Uj{9jx) and 6j = 9, where 
now 1 > > /i^^/^. One writes 

Diu~ q{x,D')u ^F + G 
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where G arises from the error term [p — Pj)uj. The bound (|2.2p is a consequence 
of the following bound (for a global £ > 0) 

(2.3) ||li||L£^L.,(|.,|<e) < f^^'O'^^P-'^^ ( Ml^lHS) + WFWl^lHS) 

+ M^6'3||(^3x2)"^m|U2(5) +/i-^6'7^||(/^^x2)2G||L2(s)) , 

and for 9 = j-i^^ 

(2.4) \\u\\Li^Ll,(\.,\<e) < A^^^"^^-'^ ( ML-ms) + \\F + GWl^l^S) ) • 

The solution u is written as a superposition of terms, each of which is product of 
Xii^i), for an interval / C [—£,£], with a functions whose wave-packet transform 
is invariant under the Hamiltonian flow of g(a:,^'). The wave-packet transform, 
which acts in the x' variables, is a simple modification of the Gaussian transform 
used by Tataru to establish Strichartz estimates for rough metrics; see also [50] , 
Precisely, set 

{T,f) ix', a - m"/' / e-<«^^'--'> {y' ~ x')) f{y') dy' . 

The base function g is taken to be of Schwartz class with g supported in a ball of 
small radius. Thus, u(x,^') = [T^w(a;i, ■)]{x\^') has the same localization in ^' as 
does u{xi,^'). 

By Lemma 4.4 of [5S] one can write 

(di - d^,q{x, a ■ d^' + d^'q{x, r) ■ d^,y{x, C) = F{x, C) + G{x, C) ■ 

By variation of parameters and the use of V2 spaces, one reduces matters to estab- 
lishing estimates for solutions invariant under the flow. The use of the V2 spaces 
from ^ 15j requires p > 2, which is implied by the conditions of Theorem ll.il 

Let Qt,s denote the Hamiltonian flow of q{x,S,'), from xi = s to xi = t. Then 
the bounds (|2.3p - (|2.4p are consequences of the following, which is the analogue of 
Theorem 7.2 of [25,. 

Theorem 2.1. Suppose that f G L^(]R^") is supported in a set of the form 
{£, ■■ Cn+i ~ M JCjl < CM , j = 2, . . . , n - 1, and ^„ sa Ofi} 

or 

{f : fn+i « Ai, ICil < c/i, j = 2, . . . ,n- 1, and |^„| < fii} 
in case 9 = fi^^^^ . 

If Wf{xi,x') = T*[f o Oo.xi] {x') , then for admissible {p, q, 7) 

Proof. The function Wf is frequency localized to ^„ « and « /i6' (respectively 
\£,n\ < when 9 sa /i^^). By duality, it suffices to show the estimate 

(2.5) \\WW*F\\l.l. < l^^'O'^^P-'^^WFh,.^,,. 
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for ^'-frequency localized F. We use t and s in place of xi and yi for ease of notation. 
Then the operator WW* applied to ^'-localized F agrees with integration against 
the kernel 

Kit, x'- s, y') = ^.t / e'(Cy-^)-(C.,*,y'-^..)5(^^ (3,' _ g^^i (y' _ ^^^^)) (^) dz dC 



where (zs^tTCs.t) — 0s.t(-z,C)- To align with the notation that x' = {x2, ■ ■ ■ ,Xn+i) 
denote the space parameters, we take C, — {(2, ■ ■ ■ , Cn+i)- Then /3e(C) is a smooth 
cutoff to the set 

{C : Cn+i « M JCjl < c/i ,i = 2, . . . , n - 1, and Cn ~ Ofi} 
(respectively |C„| < /i3 in case 9 = 

Analogous to [23 ('''•l)-(7.2)], we establish the inequalities 



(2.6) 
and 

(2.7) I K{t,x';s,y')f{y')dy 
Interpolation then yields that 
j K{t,x';s,y')f{y')dy' 



K{t,x'-s,y')f{y')dy' 



< 



< 



/i"0(l+Ai|t-s|)-'^(l+A*0'|i-s|)-^||/|Li, 



< (/i"0)'^' (1 + M|t - s\y'^^'~l\l + fiO'lt - s|)-^(i-f )||/||^,, 

y' 

In the case ^^(1 — |) < | < ^^(1 ^ f )> ^^'^ exponent in the third factor on the 
right can be replaced by ^V^(l — -) — -< 0, showing that 



K{t,x'-s,y')f{y')dy' 



<^2,^2((„-l)(i-l)-|)|^_^|-|||^||^^. 



In the case ^^^(1 —-)>-, we can ignore the last factor and obtain the bound 



K{t,x';s,y')f{y')dy 



L", 



In both cases, the Hardy-Littlewood-Sobolev inequality then establishes ()2.5p . 

The inequality (|2.6p is estimate [551 (7-1)], which follows from the fact that 
is an isometry and Qt,s is a measure-preserving diffeomorphism. Hence it suffices 
to prove (|2.7p . As in [25], we consider two cases. 

In the case ii9'^\t — s| > 1, we fix 6* < so that ^0 \t — s\ = 1, and decompose 
/3e(C) into a sum of cutoffs Pj{C), each of which is locaHzed to a cone of angle 6 
about some direction Q. The proof of [25, Theorem 5.4] yields that 

\m,x';s,y')\ < ^,-r-\l + ^l9\y'-x',^,Jy'' , 
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where x'^^^ j is the space component of 0s,t(a;, (j). For each fixed (s, t) the Xg ^ j are 
a {ij9)~^ separated set, and adding over j yields the desired bounds, since in this 
case 

= Ai"^ \t-s\-^ < n'^e {1 + n\t - si)-"^ {1 + ^i9^\t - s\)-i 

In case fj,9^\t — s| < f , we let 9 > 9 be given by 

9 = niin(^"^|i - sp^ ^ i ) . 

FoUowing the proof of [25, (7.2)], we set C" — (C2, • • • , Cn-i,Cn+i), and let Pj be a 
partition of unity in cones of angle 9 on M"^^. We then decompose 

/39(C) = ^/39(C)/3,(C"), 
3 

and let K = Kj denote the corresponding kernel decomposition. 
The arguments on page 152 of yield 

\K,{t,x'; s, y')\ < ^"0""' 9{l + ^i9\{y'- x'^^.^^h^n-il)'"' • 

The x'g f j are {^9)^^ separated in the (2,...,n — 1) variables as j varies, and 
summing over j yields 

\K{t, x'; s, y')\ < fi''9T'^ « /i"^ (1 + mI^ - sl)""^ (1 + M^'K - . □ 

3. Applications to semilinear wave equations 

As an application, we consider the following family of semilinear wave equations 
with defocusing nonlinearity 

(3.1) dfu^Au+lur'u^O {u,dtu)\t=o = if,g) u\oM = 0, 
or 

(3.2) d^u~ Au + \u\'-\^0 {u,dtu)\t=o = if,g) d.u\aM = 0, 

We will be mostly interested in the range of exponents r < 1 + ^^^2 (energy sub- 
critical) and r = 1 + -^^^ (energy critical). 

In the boundaryless case where 17 = R", the first results for the critical wave 
equation were obtained by Grillakis [TD]. He showed that when n = 3 there are 
global smooth solutions of the critical wave equation, r = 5, if the data is smooth. 
Shatah and Struwe [50] extended his theorem by showing that there are global 
solutions for data lying in the energy space x L^. They also obtained results 
for critical wave equations in higher dimensions. 

For the case of obstacles, the first results are due to Smith and Sogge [23]. They 
showed that Grillakis' theorem extends to the case where is the complement of a 
smooth, compact, convex obstacle and Dirichlet boundary conditions are imposed, 
i.e. (|3.ip for r = 5. Recently this result was extended to the case of arbitrary 
domains in O C and data in the energy space by Burq, Lebeau and Planchon 
[5]. The case of nonlinear critical Neumann- wave equations in 3-dimensions, (j3.2p . 
was subsequently handled by Burq and Planchon [5]. 
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The proofs of the resuhs for arbitrary domains in S-dimensions used two new 
ingredients. First, the estimates of Smith and Sogge [25j for spectral clusters turned 
out to be strong enough to prove certain Strichartz estimates for the linear wave 
equations with either Dirichlet or Neumann boundary conditions. Specifially, Burq, 

r — 5 

Lebeau and Planchon [5] showed that one can control the L'^Wq"' norm of the 
solution of (1.1) over [0, 1] x in terms of the energy norm of the data, assuming 
that n is compact. The other novelty was new estimates for the restriction of u 
to the boundary, specifically Proposition 3.2 in [5] and Proposition 3.1 in [6]. In 
the earlier case of convex obstacles and Dirichlet boundary conditions treated in 
|23j such estimates were not necessary since for the flux arguments that were used 
to treat the nonlinear wave equation (|3.ip . the boundary terms had a favorable 
sign. We remark that by using the results in Theorem 1.1, we can simplify the 
arguments in [5] and [6] since we now have control of the LfL^^([0, 1] x norms 
of the solution of (1.1) in terms of the energy norm of the data. If this is combined 
with the aforementioned boundary estimates in [5] and [3] one can prove the global 
existence results in these papers by using the now-standard arguments that are 
found in [23] for convex obstacles, and [20j and [26] for the case where = R^. In 
the next section we shall show how these LfL^,^ and the weaker L^L^P estimates 
can be used to show that there is scattering for (j3.ip when n — 3, r = 5 and fl is 
the compliment of a star-shaped obstacle. 

Let us conclude this section by presenting another new result. We shall show that 
the Strichartz estimates in Theorem 1.1 are strong enough to prove the following: 

Theorem 3.1. Suppose that ft cW^ is a domain with smooth compact boundary. 
IJKr <i and {f,g) G {H^{n) D L''+^{n)) x L^{n) then and ([321) have a 

unique global solution satisfying 

ueC"{[o,T];H^{n)nL''+\n))nc\[o,T];L^{n))nL^tLl{[o,T] X n) 

for every T > 0. Ifr — 3 then the same result holds provided that the {H^CiL'^) x 
norm of (/, g) is sufficiently small. 

The local existence results follow from the fact that Theorem 1 1 . 1 1 implies that if 
{df — A)v = F and v has either Dirichlet or Neumann boundary conditions then 
for < r < 1 there is a constant C so that 

(3.3) \\v\\^LUio.T)xn)<c[\\v{0, ■)\\h^ + \\dtv{0, + \\F{s, ■)hds). 

If is the complement of a bounded set, then estimate (3.3) holds with replaced 
by , as can be seen by combining the estimates for the case of compact O with 
the global Strichartz estimates on M^, and using finite propagation velocity. Using 
this estimate the theorem follows from a standard convergent iteration argument 
with u in the space 

X = c°((o, T); H\n) n u-+^{n)) n c\{q, T)-L^{n)) n i?i^((o, t) x n), 

and T being sufficiently small depending on the {H^ DL^'^^) x norm of the initial 
data {f,g) of either (3.1) or (3.2) for 1 < ?■ < 3, and T depending on the data in 
the critical case r = 3. For data of sufficiently small norm, one can obtain existence 
for T = 1 for the critical case r = 3. Together with energy conservation, the above 
yields global existence for 1 < r < 3, and global existence for small data for ?■ = 3. 
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The analog of ()3.3p when n = 3 involves L^L^P in the left. As we mentioned 
before, a stronger inequality involving L^L]^ is valid when rt = 3 by Theorem ll.il 
Any sueh corresponding improvement of (|3.3p when n = 4 would lead to a global 
existence theorem for arbitrary data for the critical case where r = 3, but, at 
present, we are unable to obtain such a result. 

4. Scattering for star-shaped obstacles in 3-dimensions 

We now consider solutions to the energy critical nonlinear wave equation in 3+1 
dimensions in a domain 17 = \ /C exterior to a compact, nontrapping obstacle /C 
with smooth boundary 

Uu{t, x) = {df - A)u{t, x) = -u^it, x), {t,x) eRxD, 

(4-1) "Lxao^O 

{Vu{t,-),dtu{t,-)) e L^{n) teR 

We restrict attention to real- valued solutions u{t, x). 

When /C is a nontrapping obstacle, the estimates above, combined with those of 
Smith and Sogge [M] (see also Burq [3], Metcalfe [T7]) imply the following estimate 
on functions w{t, x) satisfying homogeneous Dirichlet boundary conditions 

(4.2) ||u'||L£i(R;Li"(0)) + ||w||l4(M;L12(0)) 

< C(|| (V,U.(0,-),9tw(0,-))||L^(O) + |lnU'|lLi(R;L2(n))) • 

In this section, we show how these global estimates can be used to show that solu- 
tions to the nonlinear equation (|4.ip above scatter to a solution to the homogeneous 
equation 

av{t,x) = 0, {t,x)eRxn 

(4-3) HExao = 

{Vv{t,-),dtv{t,-)) e L^{n) teR. 

Let v — i^{x) denote the outward pointing unit normal vector to the boundary at 
X € die. We call the obstacle IC star-shaped with respect to the origin if ^{x) ■ x > 
for all X e dJC. Define the energy functional 

Eo{v;t)^^ ( \W.,v{t,x)\'' + \dtv{t,x)\^ dx, 

and recall that t ^ EQ{v;t) is conserved whenever is a solution to the homoge- 
neous equation (|4.3p . We show the following: 



Proposition 4.1. Suppose u solves the nonlinear problem (|4.ip and that JC is star- 
shaped with respect to the origin. Then there exists unique solutions v± to (|4.3p 
such that 

(4.4) lim Eo{u~v±;t) =0. 

t — >±oo 

Moreover, u satisfies the space-time integrability bound 

(4.5) ll'"||L5(R.iio(o)) + ||'"|U4(B;Li2(n)) < oo. 
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When /C = 0, this follows from the observations of Bahouri and Gerard [T]. 
We also remark that when /C is convex, similar results for compactly supported, 
subcritical nonlinearities were obtained by Bchatnia and Daoulatli [3]. 

Attention will be restricted to the function, as symmetric arguments will 
yield the existence of a v- asymptotic to u at — oo. As observed in [1], we actually 
have that (|4.4p follows as a consequence of (14. 5p . We first establish the existence 
of the wave operator, namely that for any solution v to (j4.3p . there exists a unique 
solution u to (|4.ip such that 

lim Eo{u — v;t) — 0. 

t — >oo 

Given (14. 2p . for any 6 > we may select T large so that \\v\\L5(^[T,oo)-L^o(n)) < S. 
Given any w{t,x) satisfying ||w||i5([7- oo);L1"(si)) ^ we have a unique solution to 
the linear problem 

Dw ^ -{v + w)^ 
lim Eo{w;t) = 

t — >oo 

as the right hand side is in L^{[T, oo); L^{il)). The estimate (|4.2[) then also ensures 
that 

I|w||l5([t,oo);L"'(0)) < C'|iw + w||i5([T,oo);Ll0(O)) < 32C(5^ 

Hence for S sufficiently small, the map w w is seen to be a contraction on the 
ball of radius S in L^{[T,oo); L^'^{fl)). The unique fixed point w can be uniquely 
extended over all of M x fi. Hence taking u — v + w shows existence of the wave 
operator. 

To see that the wave operator is surjective, we need a decay estimate which 
establishes that the nonlinear effects of the solution map for (|4.ip diminish as time 
evolves. 

Lemma 4.2. Let tC be star-shaped with respect to the origin. Ifu{t,x) solves (|4.ip . 
then the following decay estimate holds 

lim I / \u{t,x)\^ dx = 0. 

t-*oo 6 

When JC — 0, this is due to Bahouri and Shatah [2 . The proof below is es- 
sentially theirs, with slight modifications made to handle the boundary conditions. 
However, for the sake of completeness, we replicate the full proof below. We re- 
mark that the approach has its roots in arguments of Morawetz , and is related 
to other wrks regarding the decay of local energy for linear solutions in domains 
exterior to a star-shaped obstacle. 

To see that this implies the proposition, observe that given any e > 0, there 
exists T sufficiently large such that 

sup||w(i, OIIls < s. 

t>T 
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Hence for any S > T we obtain the following for any solution u to (|4.ip 

\\u\\L'^i[T,S]:L^O{n)) + \\u\\L''{[T,S];L^^{n)) < C {E + \\u^\\L^[T,S]:L'^{n))) 

< CE + Ce\\u\\Li^iT,s];L^''{ny) 

where E denotes the conserved quantity 

E ^ E{t) ^ J^^^\Vuit,x)\-' + ^\dtuit,x)\^ + ^\u{t,x)\'' dx. 

A continuity argument now yields ||M|lL5([T,oo);Lio(f2)) + l|u||Li([T,cx3);Li2(n)) < '^GE 
and by a time reflection argument, (|4.5p follows. However, this implies that the 
linear problem 

Uw^-u^ \miEQ{w]t)^Q 

t — ^oo 

admits a solution, showing that the wave operator is indeed surjective as v = u — w 
is the desired solution to (|4.3p . 

Proof of Lemma \4-'^ By a limiting argument it suffices to consider smooth, classical 
solutions u which decay at infinity. We must show that for for any Eq > 0, there 
exists To such that whenever t > Tq, 

'u{t, x)\ dx < Eq- 



6 

Consider the stress energy tensor associated with u (see Tao [3S], p. 149) 

T°^ = -dtud^^u l<j<3 

T'" = d^^ud^.u ~ ^(|Vup - (dtu)^ + iyS) 1 < J, fc < 3. 

It can be checked that the divergence free property holds 

dtT°° + d^^T^^ ^ (9tT"^ + d.^.P'' = 

with the summation convention in effect. Taking the first of these identities and 
applying the divergence theorem to a region {0 < t < T,\x\ > R + t} (with R > 
large enough so that K, C 5^(0)) we have 

(4.6) / l\dtuiT,x)\^ + l\Vu{T,x)f + l\uiT,x)fdx + ^&u^iO,T) 

J\x\>R+T Z Z b V2 

< / i|9,^(0,a;)p + i|Vu(0,x)|2 + i|u(0,x)|6dx 
J\x\>R z z b 



where 







/ 1 
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M§ := {a <t <b,\x\ ^ R + t} 

Since the solution has finite energy, we may select R large so that the right hand 
side of (|4.6p is less than |^ (and again JC C Bii{0)). By time translation, 1 1-^ t + R, 
it will suffice to show the existence of Tq such that whenever t > Tq we have 



I I \u{t,x)\ dx<^ 
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(the additional smallness in the right hand side of ()4.6p will be used later in the 
proof). 

We now define the following vector field X — {X'^ , , X'^ , X^) by contracting 
the stress-energy tensor with the null vector field tdt — x-Vx and adding a correction 
term 

X" = tT"" - XkT"'' + udtu 
X^ = tT^° - XkT^^ ~ ud^^u l<j<3. 
The space-time divergence of X satisfies 



div(X) = - 



1 



3 

We now apply the divergence theorem over the truncated cone K^l={x&n: 
\x\<t,Ti<t<T2} 

XUx~ [ XUx- [ lx"-y^xAda 

D{T,) Jd(T,) J Mil \ 1^1 I""' ) 

+ [ ^-^dxdt-f 1^ ■ {X\X^,X^) da 

= I + II + III + IV + v 

where da denotes Lebesgue measure on the corresponding surface and D{Ti) = {x d 
n : \x\ < rj, and M^^ = {\x\ = t,Ti < t < T2}. The star-shaped assumption 
is crucial in controlling the last term V. Indeed, consider the restriction of the 
integrand in V to the boundary d^l{= dJC) and observe that the Dirichlet boundary 
condition gives 



v{X\X\X^) = - ^iXk[dx^udx,u-^\\/u\ 

l<j,k<3 



I 



2 

2 



= -{iy Vu) (x • Vm) + 2 ■ a;)|V'u| 

We have that Vu is normal to dfl and hence |Vup = {v ■ Vuf'. Treating 
vector, we can project it on to the subspace orthogonal to u obtaining 

= Vu ■ [x ~ {v ■ x)v) — X ■ Vu — (y ■ x){v ■ Vu). 

This now gives 

V ■ {X\X^, X^) = -\{v- x) {v ■ Vuf < 
and since IV > is clear, 

> / + // + ///. 

We now impose polar coordinates (r, w) G M x §^ on the third term, writing 
/// = — ^ / {r{dtu + druf + u{dtu + dru)) da 
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where dr — jfy ' ^ denotes the radial derivative. Next parameterize Af^^ by (r, uj) — > 
(r, roj) and set v{y) — u{\y\,y) (or v(ruj) = u{r,ruj) in polar coordinates) so that 
we may write concisely 

III ^- / r fa^w + -) drduj+ / -9^ (r^^) rf^ 

- / / r(drV + -Y drdLj + - [ T^v^ (T2UJ) dw - - [ T?v^{Tiu;) duj 
To handle the first term /, first observe that in polar coordinates 

Since /C is star-shaped we may parameterize dfl by (r, cj) = (\E'(a;),a;) where $ is a 
real valued function on This allows us to write 



Jd(T2) ^ \ 



1=1 ^1 (dtuf +{dru+ ^u] + ;^|V„up + iwM + r(dr + ^u] dtu dx 



(4.7) 

1 

~ 2 



§2 



T2 

T2dr{ru'^) dr dto 



Integrating by parts in the last term yields cancellation with one of the terms in /// 
as the boundary condition gives — i Jga /^^^-j T2dr{rv?) dr duj — — ^ T2V^(T2Lu) dio. 
Similarly, 



/ ^({dtuf+(drU+-l 

Join) ^ \ \ r 



II = - I ^1 (dtu)'' +{dru+-u] + ^\V^u\^ + ^uA + rfdr + ^u] dtU dx 



+ \ I TyiTiiu) duo 

In order to control remaining term in / we need to observe the following Hardy 
inequality, which holds in the exterior domain 

(4.8) / ^ < 4 / |Vup dx. 

Jn \A Jn 

To see this, we assume u is real-valued and denote the integral on left hand side as 
J. Converting to polar coordinates 

J = / dr duj — / ru{ru!)'^ du) — / 2u{dru)r dr dui. 



S2 J-i>(uj) JS^ 



The first term on the right is nonpositive (provided u exhibits sufficient decay at 
infinity) and Cauchy-Schwartz on the second term gives 



J < 2VJ ( [ I la^upr^ dr duj 



The inequality (|4.8p now follows 
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We now observe that the first integral in (|4.7p is bounded below by T2 /oj-j-^) dx. 
Setting Ta = T > and Ti eT (0 < e < 1) and using the Hardy inequality (HiHD 
to control the first integral in // now yields 

T [ dx< CeTE + [ [ T (drV + -Y duj dr. 

Jd{T) 6 JeTJS2 V rJ 

Here E is the conserved quantity E = E{t) = J^^ T^^{t^ x) dx. We can now divide 
both sides of this inequality by T and choose e sufficiently small so that CeE < Eq/A, 
leaving us to control the integral involving v. However, by the proof of the Hardy 
inequality above we have 

/ f (drV+ -Y duj dr < 10 [ [ {drv f r^ duj dr < 10 R\ix{eT, 00) < -, 
JeT V rJ Js^ 4 

provided T is large enough so that eT > R. □ 
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